% 来自文献DOI: 10.5120/11756-4989的模型
function out = twta(input)

inputPha = angle(input);
inputMag = abs(input);

aa = 2;
ba = 1;
ap = 4;
bp = 9;

MagDist = aa * inputMag ./ (1 + ba * inputMag.^2);
PhaDist = ap * inputMag.^2 ./ (1 + bp * inputMag.^2);

out = MagDist .* exp(1i * (inputPha + PhaDist));
end